Dipolar bosons in a planar array of one-dimensional tubes 
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We investigate bosonic atoms or molecules interacting via dipolar interactions in a planar array of 
one-dimensional tubes. We consider the situation in which the dipoles are oriented perpendicular to 
the tubes by an external field. We find various quantum phases reaching from a 'sliding Luttinger 
liquid' phase to a two-dimensional charge density wave ordered phase. Two difi'erent kinds of charge 
density wave order occur: a stripe phase in which the bosons in different tubes are aligned and a 
checkerboard phase. We further point out how to distinguish the occurring phases experimentally. 
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Recent experiments in ultracold atoms have achieved 
the realization of various quantum phases. These phases 
reach from the superfluid and Mott-insulator of bosons 
in optical lattices to a BCS phase and a Bose-Einstcin 
condensate of molecules in fermionic gases [H. In addi- 
tion, it has been possible to engineer different trapping 
geometries which enabled the study of the remarkable 
physics of one-dimensional quantum systems 0, H, 01 ■ In 
such systems, known as Luttinger liquids the interac- 
tions play a major role and lead to properties quite differ- 
ent from their higher-dimensional counterparts. This has 
been nicely demonstrated investigating the crossover be- 
tween an array of decoupled one-dimensional tubes and a 
quasi-three-dimensional system by lowering the potential 
barrier between the tubes i^. Since the interactions in 
most cold atomic gases have a short range of the order of 
a few nm and the distance between the tubes is typically 
500 nm, the only coupling between the tubes is gener- 
ally provided by the hopping of particles between tubes. 
For bosons this leads to a dimensional crossover between 
the peculiar one-dimensional phases and a quasi-three- 
dimensional superfluid 

However, other fascinating phases may occur when 
one-dimensional tubes are directly coupled by interac- 
tions. These include in particular the so called 'slid- 
ing Luttinger hquid' (SLL) 0, H, that was studied 
in connection with high-Tc superconductors and stripe 
physics for fermionic systems. In the SLL phase the 
typical properties of one-dimensional systems, namely al- 
gebraically decaying correlations, survive despite the cou- 
pling. Thus the Fermi liquid phase which usually occurs 
in more than one dimension is suppressed. Experimen- 
tally the SLL has not been observed yet. 

In cold atomic gases these phases could not be ex- 
plored so far due to the lack of interactions extend- 
ing over the range of the inter-tube distances. How- 
ever the experimental progress in the realization of 
quantum degenerate atomic and molecular gases with 
dominating dipole-dipole interactions llj have shown 
the potential to bridge that gap. Particularly promis- 
ing are polar molecules with large electric dipole mo- 
ments [H El, [ll, [S [3 , since the strength of the dipolar 
interaction is considerable over the range of tube spac- 



ings [Trj . 

In this manuscript we investigate the possible quantum 
phases of a dipolar bosonic gas in a planar array of one- 
dimensional homogeneous tubes (Fig. [1]) . This situation 
can be achieved experimentally by using a strong two- 
dimensional optical lattice. If sufhciently strong, this lat- 
tice suppresses the hopping of particles between different 
tubes. Even in this situation the dipole-dipole interaction 
couples the tubes. The dipole-dipole interaction is attrac- 
tive or repulsive depending on the relative orientation of 
the dipoles. This anisotropy causes interesting phenom- 
ena such as an instability towards collapse [l^, |T3| ■ Here 
we focus on the situation where the dipoles are aligned 
by an additional external field. The orientation is cho- 
sen perpendicular to the direction of the tubes, since this 
configuration is the most stable [23|. The orientation 
with respect to the plane of the array is varied. 




FIG. 1: A planar array of one-dimensional tubes. The dipoles 
point perpendicular to the tube direction. In the inset the 
angle a of the orientation of the dipoles with respect to the 
plane of the array is shown. 

Assuming vanishing hopping between the tubes we 
find four different regimes: a SLL state with (i) domi- 
nating superfluid correlations or (ii) dominating charge 
density wave (CDW) correlations along the tube forms, 
(iii) checkerboard or stripe CDW ordering develops 
(cf. Fig. [3]), and (iv) an instability towards collapse occurs 
[28| . The finding of a SLL phase is novel in a bosonic sys- 
tem. The SLL has power-law correlation functions along 
the tubes which can be dominating superfluid correla- 
tions (i) or dominating CDW correlations (ii). The two 
regimes are connected by a crossover. We discuss how 



2 



cold quantum gases could provide the opportunity to ob- 
serve the SLL phase for the first time experimentally. 

Interacting bosons in an array of one-dimensional 
tubes (cf. Fig. [T]) can be described by the Hamiltonian 

j 

+ E Jdzdz'Vd{{j-j')a,0,z~z')p,{z)p,,{z').{l) 

Here is the bosonic anihiliation operator in tube j and 

Pj = ^j^j is the density operator. We use h = 1. The 
tube distance is a, M is the mass of the particles, and g 
is the strength of the ^-interaction between the particles. 

The last term in Hamiltonian ([T]) stems from the 
dipole-dipole interaction. Its amplitude is Vd{r) = 
Vb(d^ — 3(df)^)/r'^, where d = (cos a, sin a, 0) denotes 
the direction of the dipole moment and f is the unit 
vector. The interaction strength Vq is given by Vq = 
cPp,o/{4:TT) for magnetic and Vq = d^/(47reo) for electric 
dipoles. Here d is the strength of the dipole moment, and 
/j,o (eo) is the vacuum permeability (permittivity). 

The low-energy properties of the system can be de- 
scribed using the bosonization approach 0, 0|. Two 
bosonic fields, 4'j and 0j, are introduced to describe the 
modes of the system, where (j>j and 9j are related to the 
amplitude and phase of the operator '^j, respectively. 
The Hamiltonian in the absence of the dipolar interac- 
tion becomes 

^" = |^E/d^ (Kive,{z)r + ^{y^,{z))A,{2) 

where the parameters u, the velocity, and K, the Lut- 
tinger parameter, depend on the microscopic details of 
the underlying system For weak interaction the re- 
lations K w '^^J^Tg ^^"d " ^ hold. Here po is the 
average density of bosons. 

The dipolar interaction can be expressed in the 
bosonization language using pj{z) = po — (l/7r)V0j(z) + 
PoEp/oe'^''^'''"'^~'^'^''^^ where p is an integer. The 
smooth part of the density ~ V(j)j yields a non-local 
quadratic term in the Hamiltonian whereas the first har- 
monic \p\ = 1 introduces backscattering terms. Neglect- 
ing higher harmonics, \p\ > 1, the contribution to the 
Hamiltonian of the dipolar interaction thus becomes 

Hd = — 2 E / d^;dz' Vd(na, 0, z-z') X 

X (v^,(z)V0,+„(z') -I- O^^iz, z')) , (3) 

where Oj„(z, z') oc cos{2(j)j{z) — 2(j)j^n{z')) are the back- 
ward scattering operators. 



Assuming that the average particle spacing is large 
compared to the size of the dipoles, the detailed struc- 
ture of the interaction at short distances can be ne- 
glected [20!, '2l!|, and any possible further contribution 
to the interaction on a length scale below a cutoff 1/po 
can be absorbed into the values of u and K. By this 
definition K can now range from very large values for 
a weakly interacting gas without dipolar interactions to 
values smaller than one in t he p resence of finite dipolar 
interactions as shown in Ref. [221 for a purely dipolar one- 
dimensional system. Reaching values of K smaller than 
one in a bosonic system is a characteristic of non-local in- 
teractions Q . In the following we use the dimensionless 
ratio of the dipolar interactions and the s-wave scattering 
given by 7 = 4Vbi^/(a^u). 

In order to determine the quantum phases occur- 
ring in the system, we investigate the relevance of 
the superfluid and CDW operators along and perpen- 
dicular to the tubes. To find the dominant corre- 
lations we first consider the effects of the quadratic 
part of Hq + Hd by comparing the relevance of the 
density-density correlations {pj{z)pj{0)) oc (l/z)^^'-^*.ii 

and the single particle correlations (^|(z)*l'j(0)) oc 
(l/z)"-^'''ii/(^-^). The algebraic decay of both correla- 
tions is typical for a Luttinger liquid Q. We find 
that the exponents depend on the strength of the 
dipolar interaction and the orientation of the dipoles 
through G^,||(7,a) = l/{2n)j:,J^ + lV{q^,0)r'/^ 
and Ge,||(7,a) = 1/ {2n)j:^Jl + jV{q^,0)y/\ Here 
V{q±,q\\ = 0) = 2(a/9o)2-4cos(2a)Re[Li2(e'''«^)], where 

^(9-L> ^ll) ^d(9±,9||)/(47rVo) denotes the dimensionless 
interaction strength, Li2 the polylogarithm, and Re the 
real part. q± and g|| are the components of the momen- 
tum perpendicular and parallel to the tubes, respectively. 
To perform the calculations a discrete Fourier transform 
is taken in the a;-direction perpendicular to the tubes, 
and a continuous Fourier transform along the tubes. We 
further approximate the momentum dependence of the 
dipolar interaction along the tubes by the low-momentum 
value at q\\ = which yields the main contribution to the 
long-distance behavior of the correlation functions. This 
corresponds to using an effective interaction that is local 
along the tube and whose strength is given by the dipolar 
interaction integrated over the tube [23 ]. 

Comparing the obtained exponents, the CDW correla- 
tions in the tube dominate at long distances if 

2XG0,|| - ^Ggji < 0. 

For decoupled tubes the functions reduce to || = 1 
and Gen ~ ^^^^ case one recovers the crossover 

for a single tube [5, 22] between a region in which the 
superfluid correlations are dominant for i^T > 1/2 and a 
region in which the CDW correlations are dominant for 
K <l/2 (cf. 7 = line in Fig. [2]). 

For coupled tubes the inter-tube operators can be- 
come relevant and change the nature of the occurring 
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FIG. 2: (Color online) Different quantum phases occuring for the directions (a) a = and (b) a = tt/2. Experimentally K can 
be changed varying the s-wave scattering length of the particles. Further tuning the dipolar interaction would correspond to 
changing both 7 and K. The subscript 'z' denotes the dominant correlations along the tube. 



quantum phases. If the pre-factors of the operators 
„ are small, their relevance can be determined by 
looking at their scaling dimension. The scaling dimen- 
sion of the operators is 2KG^^j_^n using G^^j^^n = 

1/(2^) (1 + lV{q±,0))-^/^ (1 - cos{naq±)). There- 
fore a backscattering operator between tubes at distance 
na becomes relevant and can induce charge ordering, if 

2KG^,±^n - 2 < 0. 

In Fig. [5] the phase diagram depending on 7 and K is 
shown for two different orientations of the dipoles [30| . 
Experimentally K can be varied by adjusting the short- 
range interaction between the particles. Varying the 
dipolar interaction strength changes both, 7 and K . In 
Fig. [2la) the dipoles are pointing along the direction of 
the array {a = 0) leading to an attractive inter-tube in- 
teraction at short distances. By contrast in Fig. [H^b) the 
dipoles are oriented perpendicular to the plane of the ar- 
ray (a = 7r/2) resulting in a purely repulsive inter-tube 
interaction. The phase boundaries are shown for the case 
where the spacing between tubes equals the average par- 
ticle distance, i.e., apo = 1 [sij. 

For 7 <C 1 the behavior for the different orientations of 
the dipoles is very similar: For K > 1 a SLL with dom- 
inant superfluid correlations (regime (i)) occurs. In con- 
trast to an array of independent tubes, the presence of the 
dipolar interaction leads to a coupling of the densities in 
different tubes on long length scales. In the bosonization 
language the coupling is given by the forward scattering 
terms, i.e., in Hamiltonian ([3]) the terms which contain 

V0j(z)V0,v(2') with j ^ j'. 

The CDW correlations along the tubes become rele- 
vant for values K < 1/2. However, the CDW correlations 
perpendicular to the tubes become relevant already for 
larger values K < 1. In the regime 1/2 < K < 1 (regime 
(iiia)), a recalculation of the intra-tube correlations tak- 
ing into account the ordering in the perpendicular direc- 
tion then leads to a CDW ordered phase in all directions. 
In this case, the order in the tube will be weaker than 



the order in the direction perpendicular to the tubes. In 
contrast, for K < 1/2, the order along the tubes is ap- 
proximately as strong as the order perpendicular to the 
tubes, since both the CDW ordering along and perpen- 
dicular to the tubes are relevant (regime (iiib)). 

Above a critical value of 7 an instability towards col- 
lapse (regime (iv)) occurs. The instability occurs if there 
exists a q± for which 1 -I- ^V{q±) < 0. The physics of 
the instability is different for the shown cases a = 
and a = tt/2. For a = the attractive inter-tube in- 
teraction overcomes the repulsive intra-tube interation 
and causes the collapse of the particles inside the tubes 
and an alignement of these perpendicular to the tubes. 
However, the instability for a ~ tt/2 stems from a strong 
repulsive inter-tube interaction which dominates over the 
intra-tube interaction. This can occur, if the tubes are 
very close and the contact interaction is attractive. Here 
the particles collapse inside a tube in order to avoid the 
interaction with the particles in the neighbouring tubes. 

For intermediate values of 7 the behavior depends on 
the orientation of the dipoles. For a = the CDW or- 
der can already be reached for values oi K > 1, i.e., the 
dipolar interaction enhances the tendency of the system 
to order. By contrast for a = n/2, a larger value of 7 
can destabilize the CDW ordering. Thereby a transition 
between a SLL and a CDW ordered phase seems to be 
possible simply by varying the orientation of the dipoles. 
Further for a = 7r/2 a small region of SLL with dominat- 
ing CDW order (regime (ii)) can be seen in Fig. ^h) at 
large 7. Whether this survives for realistic experimental 
parameters is an open question. 

In phase (iii), the form of the CDW order which oc- 
curs depends on the direction of the dipoles (Fig. [3]). In 
particular if the dipoles lie in the plane of the tubes, 
the interaction between tubes is attractive at short dis- 
tances and therefore the CDW order of different tubes 
is aligned [1^. By contrast, if the dipoles are perpen- 
dicular to the plane, the interaction between tubes is 
repulsive and a checkerboard ordering is expected. To 
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FIG. 3: A sketch of the transition between a stripe and 
checkerboard CDW order depending on the orientation of the 
dipoles is shown. 



determine which form of CDW order the system takes 
for a given angle a we perform a minimization of the 
energy of the system assuming CDW order afong the 
tubes, i.e., we take the density in tube j to be of the form 
pj{z) = po (l + cos{2npoz — 4>j))- Here 0j describes the 
average phase in tube j. With this ansatz (and using 
5z = z — z') the expression to be minimized reduces to 



{(f>j - (f)j,)). 



^ j d5zVd{{j-j')a,0,5z)coa{2TTpoSz -2 

3=^3' 



As expected, we find a transition between the stripe 
order for small a and the checkerboard order for a 
close to 7r/2 (see Fig. Assuming only couphng be- 
tween nearest-neighbor tubes the transition takes place 
at cos^(ac) — Ki{2TTapo) / {2napoK2{2Trapo)), where Ki 
are modified Bessel functions. At the transition point the 
tubes experience only a weak coupling to other tubes, 
and the correlations in the tube are superfluid or CDW 



dominated depending on the parameter regime. Taking 
the full dipolar interaction into account we determine 
the transition point numerically [s^ . The result can 
hardly be distinguished from the transition point found 
for nearest-neighbor coupling only. 

The precise setup to observe the quantum phases ex- 
perimentally depends on the realization of the dipolar 
particles. However, here we describe some of the ba- 
sic characteristics of the phases which could be detected. 
The coupling of the tubes in the SLL phase can distin- 
guish it from the Luttinger phase of decoupled tubes. 
This coupling could, e.g., be detected by exciting the 
dipole mode for part of the tubes and detecting the in- 
duced center of mass momentum of the remaining tubes. 
Further the frequency of the dipolar mode can give infor- 
mation on the state inside the tubes (2j| . The stripe and 



checkerboard orders show characteristic density-density 
correlations which could be detected measuring the noise- 
correlation spectrum of time-of- flight images [1^. Be- 
fore performing the time-of-flight measurement molecules 
could be dissociated while freezing their position by an 
additional strong optical lattice. 
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